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Abstract We calculate NLTE models of stellar winds 
of hot compact stars (central stars of planetary nebulae 
and subdwarf stars) . The studied range of subdwarf pa- 
rameters is selected to cover a large part of these stars. 
The models predict the wind hydrodynamical structure 
and provide mass-loss rates for different abundances. 
Our models show that CNO elements are important 
drivers of subdwarf winds, especially for low-luminosity 
stars. We study the effect of X-rays and instabilities on 
these winds. Due to the line-driven wind instability, a 
significant part of the wind could be very hot. 

Keywords stars: winds, outflows; stars: mass-loss; 
stars: early-type; stars: horizontal-branch; planetary 
nebulae: general 



1 Subdwarf stars and the stellar wind 

Despite being low-luminous, the horizontal branch stars 
may still have stellar winds due to their large effec- 
tive temperatures and low stellar radii. The winds of 
these stars is thus siniilar to those of luminous OB stars 
(|Vink fc Cassisill2002tlUnglaubll2008l) . which arc driven 



e.g., C, N, 0, or Fe (see 


Owockill2004; 


Krticka & Kubat 


2007; 


Puis et al. 


2008 


for a review: or HamannI 2010). 



The main wind parameters are the mass-loss rate 
M and the terminal velocity Voo- The mass- loss 
rate d epends on the s tellar luminosity L as M 



L^/" ( Puis et all 120081 ). where a' « 0.6 for O stars 
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( Krticka fc Kubatll2009al ). This scaling implies that the 
mass-loss rate increases with luminosity. Because the 
wind is driven by metal lines, the mass-loss rate de- 
pends also on the mass fraction of heavier elements, 
Z. For O star s the scaling M ~ is roughly valid 

( Krtickal 20061 ). In contrast, the wind terminal velocity 
is proportional to the surface escape speed, Voo oc Vcsc- 
Here we compl ement the stud i es of hot subdwarf 
winds provided by Vink fc Cassisi (2002) and Unglaubl 
(I2OO8I) . 



2 Wind models 

The stellar radiative flux (at the inner wind boundary) 
is taken from static, spherically symmetr ic NLTE mode l 
atmospheres calculated with the code of iKubatI (I2OO3I) . 

For calculati on of wind model s we ap ply the NLTE 
wind code of iKrticka fc KubatI boo4 . The wind 



is modelled assuming a spherically symmetric, sta- 
tionary flow. Ionization and excitation are derived 
from the solution of the statistical equilibrium equa- 
tions. Ionic models adopted from the TLUSTY 
code (jLanz fc Hubenvl 120071 ) are based on data froni 



the Opacity and Ir on Projects (jSeaton et al.l 11992 



Hummer et al.lll993() 



The derived excitation and ionization are used to cal- 
culate the radiative force and radiative heating. These 
terms are inserted in the corresponding hydrodynamical 
equations (continuity, momentum, and energy). These 
equations may be solved for each component of the flow 
separately. 

From the obtained wind structure we calculate the 
mass-loss rate and the terminal velocity. These values 
can be compared with parameters derived from obser- 
vations. 
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Fig. 1 The position of tlie studied stars in tiie HR diagram, 
compared to luminous O stars 

3 Stellar winds of central stars of planetary 
nebulae 



Up to now, the wind models of iKrticka fc KubatI (j2004^ 
were applied only to stellar winds of luminous O stars. 
Because there is a large gap between the parameters 
of luminous O stars and O subdwarfs (sec Fig. [T]), wc 
first decided to calculate wind models of central stars 
of planetary nebulae (CSPN). This enables us to test 
the applicability of our wind models for stars with low 
luminosity. 



For this test we chose the dataset of Pauldrach et al 



(2Qq3), who derived the stellar and wind parameters of 
selected CPSNs from observations. 

A comparison between wind mass-l oss rates and ter- 



minal velocities derived by us and by iPauldrach et al 



(200J, see Fig. [5]) showed good agreement. 

The contribution of individual elements to the ra- 
diative force (Fig. [3]) shows that iron contributes to the 
radiative force particularly close to the star. Lighter 
elements are important in the outer wind regions. This 
is similar to stellar winds of luminous O stars. We 
conclude that stellar wind of CSPN resembles that of 
luminous O stars. 



4 Subdvirarf stars 

Turning now our attention to subdwarf stars, we se- 
lected a grid of effective temperatures Toff and surface 
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Fig. 2 Comparison of predicted mass-loss rates (upper 
panel) a nd terminal veloc i ties (l ower panel) with results de- 
rived by IPauldrach et all (|2004l ) for CSPNe. 
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Fig. 3 The radial dependence of the contribution from 
individual elements to the radiative force in the wind of 
CSPN IC 4593 in comparison to the wind of a Cam 
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gravity acceleration g that well covers the area of sub- 
dwarf parameters (see Fig. S]). We assumed a fixed 
stellar mass M = 0.5 M©, different values of metallic- 
ity 0.3 Zq. IZq, SZq, 10 Zq, and helium-to-hydrogen 
ratios iV(Hc)/A^(H) of 0.085, and 100 (by number). 

The low luminosity of subdwarf stars implies also 
low mass-loss rates (see Fig. [5]). Contrary to the wind 
of CSPNc. the radiative force is dominated by lines of 
lighter elements everywhere, similarly to the wind of 
early B main-sequence stars (see Fig. |6l). The wind 
terminal velocity Voo is of the order of the surface escape 
speed Vcsc (Fig. [71). 
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5 Contribution of different elements 



Fig. 4 Stellar parameters of our grid of subdwarf models 



To understand how different elements contribute to the 
radiative force in different situations, let us mention 
that the radiative force due to a particular line is given 
by the line optical depth and the radiative flux. The 
line optical depth in Sobolev approximation is 



dVr 

dr 



(1) ^ 



where fij is the oscillator strength, Ui is the number 
density of atoms in a given state, and ^ is the velocity 
gradient (jPuls et al.ll2008( ). The radiative force due to 
optically thick lines (rs > 1) does not depend on level 
populations n^, whereas the radiative force in optically 
thin lines (rg < 1) depends on level populations. 

To describe the contribution of individual lines to 
the radiative force, let us divide them into strong ones 
(typically resonance lines of dominant ionization stages 
of C, N, O), interr nediate ones (strongest iron lines), 
and weak ones (cf. IVink fc de Koteill2005l ). For high 
wind density both strong and intermediate lines are 
optically thick, consequently each such line contributes 
by a similar amount (given by the radiative flux) to 
the radiative force. However, because there are much 
more iron lines than lines of CNO elements, iron lines 
dominate at large densities. On the other hand, at low 
densities strong lines may be still optically thick (and 
consequently contribute to the radiative force), whereas 
iron lines become optically thin. As the radiative force 
by optically thin lines depends on the abundance of its 
element, which is relatively low in the case of iron, the 
contribution of iron to th e radiative force is small in a 
low-density environment ( Puis et al. 200(i Vink et al. 



200lUKrticka et al.ll2006[ ). 

The stellar winds of CSPNe or O stars have rela- 
tively high density close to the star and low density in 
the outer parts. Consequently, the iron lines dominate 
the radiative driving for small radii, whereas the lighter 
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Fig. 5 The dependence of the wind radial velocity on ra- 
dius in a selected model 




Fig. 6 The radial dependence of the contribution from 
individual elements to the radiative force in the winds of a 
subdwarf star with Toft — 35 000 K, logg = 5.5, Z = Zq, 
compared to a Ori E 
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Fig. 7 The ratio of the wind terminal velocity and the 
surface escape speed for Z = Zq 



elements arc important at large radii (see Fig. [3]). On 
the other hand, the stellar wind of subdwarf stars or 
main sequence B stars has relatively low density every- 
where. Consequently, lines of light elements dominate 
at all radii (see Fig. 15]). 



6 Multicomponent wind structure 

There are more similarities between the stellar winds of 
hot subdwarf stars and of early main-sequence B stars. 
One of them is connected with inefficient transfer of mo- 
mentum from the radiatively driven ions to the passive 
ions of hydrogen and helium. 

In line driven winds, momentum from radiation is 
transfered to heavier ions, whereas hydrogen and he- 
lium practically do not receive any momentum from 
the radiation field. Because heavier elements constitute 
only a very small part of the stellar wind (of about 1% 
in mass), momentum shall be transferred from these 
heavier ions to the passive ones. As the stellar wind 
of hot stars is ionized, Coulomb collisions provide the 
most efficient way of such transfer of momentum. 

For stars with relatively high wind-density, the trans- 
fer of momentum between individual heavier ions (de- 
noted by i) and the hydrogen-helium component (p) 
is efficient, because the velocity difference between 
these components is relatively small. However, for 
low density stellar winds the momentum transfer be- 
comes inefficient and two effects may emerge: fric- 
tional heating and decouplin g of wind components 



(e.g.. iKrticka fc KubatI 120011: lOwocki fc Pulsl 12002 



Votruba et al.ll2007l : lUnglaubl 120081) . 

The importance of these effects can be assessed from 
the value of the non-dimensional velocity difference 
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Fig. 8 Frictional heating in models of the star with 
Teft = 55 000 K, logg = 6.1, and Z = 3Zq. Top: the 
non-dimensional velocity difference Eq. ((2]). Bottom: Com- 
parison of wind temperatures on one-component and four- 
component models 



where Vi, and Vp are the velocities of the components. 



2k {iTLiTp + TOpTi) 



(3) 



is the mean thermal speed, where mj, and mp arc the 
ionic masses, and Tj. Tp the corresponding tempera- 
tures. 

For small velocity differences, Xip < 0.1, the flow 
is well coupled and multicomponent effects are negligi- 
ble. For larger velocity differences, Xip > 0.1, frictional 
heating becomes important (Fig. [5]), and for Xip > 1 
the wind components may decouple. After decoupling, 
hydrogen and helium leave the star if the decoupling 
occurs at velocities larger than the escape speed, or 
fall back on the stellar surface if the decoupling oc- 
curs at velocities lower than the escape speed (e.g. 
Votruba et ai1l2007i [2OI0I ). 

The parameter domains, in which multicomponent 
effects are important are given in Fig. |9l From our 
calculations wc conclude (in agreement with Unglaubl 
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hot. On the other hand, for stars with weaker winds 
this ratio could be relatively large, indicating that once 
a shock occurs in the wind, the material has too low 
density to cool down significantly radiatively. Conse- 
quently, for low-luminosity stars a significant portion 
of the wind can be very hot, and it can be difficult to 
observe such winds. 

For low-luminosity O stars this can explain the 
so-called "weak wind problem" (see iKrticka fc Kubat 



2009al) . The same effect could be also present in subd- 
warf stars fFig. [TU)) . 



Fig. 9 Multicomponent effects in the Toft — log g diagram 
-,3 



10" 
10^ 
10 
1 

0.1 

10-2 
10"^ 
10-4 



1 1 1 1 — 


-1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 

r3f, = 35 000 K, log 9=5.5. 








.f'^"" Tgfl = 55 000 K, log 9=5.5 ■ 












subdwarf star 




weak wind problem star 




normal star 



4 6 

r/R,. 



10 



8 Similarity to the first stars 

Stellar winds of hot subdwarf stars are also similar to 
the stellar wind of the first stars in the Universe. First 
stars formed from the pristine gas that emerged from 
the primordial nucleosynthesis. Consequently, these 
stars were pure hvdrogen- helium, and were very mas- 



IGOMq, e.g.. lOmukai fc Fallal [2OO3I ) and 



sive (M 

hot. Hydrogen and helium lines ar e not able to ac- 
celerat e the wind from these stars (|Krticka &: Kubat 
2009bl) . however, during the core He-burning phase 
the CNO elements (pri mary nitrogen ) eme rge at the 



surface due to mixing ( Mevnet et al. 20061) . Conse 



quently, similarly to subdw arfs stars, evolved firs t stars 
had CNO-driven wind s (jCrafener &: Hamann 
Krticka fc Kubatll2009bl) . 



2008 



Fig. 10 Dependence of the ratio of cooling and hydrody- 
namical time scales on radius for individual stars 



20081 ) that coupled subdwarf winds with 
M < 10"^^ Mq yr^i do not exist. 



7 Inefficient shock cooling 

The last similarity between the stellar winds of B stars 
and subdwarf stars mentioned here is connected with 
inefficient shock cooling. Hot star winds are known to 
be subject of the instability of radiative driving. As a 
result of this, a small part of the wind ma terial of lumi 



nous O stars is very hot and emits X-rays (jOwocki et al, 
1988t iFeldmeier et al.|[f997l ). 



However, for low-luminosity stars (e.g., the subd- 
warf stars) the ho t part of the stella r wind may be sig- 
nifican tly larger ( Cohen et al. 20081 : Krticka fc Kubat 
2009a[) . This can be seen from the ratio of cooling and 
hydrodynamical time scales plotted in Fig. [TOl For O 
stars of normal luminosity this ratio is very small, in- 
dicating that only a very small part of the wind is very 



9 Importance of the stellar winds 

Although stellar winds of hot subdwarf stars are rather 
weak, they can still have some observable effects. 

Weak stellar winds were invoked for the explanation 
why some subdwarf stars pulsate whereas othe rs do not 
(e.g., Fontaine et ahlboOG ). In agreement with Unglaubl 



(|2008l) , our results do not support this explanation, be- 
cause weak stellar winds are expected to be decoupled. 

On the other hand, weak stellar winds in binaries 
consisting of a subdwarf star and a compact object 
(neutron star or a black hole) may lead to production 
of X-rays. Last but not least, interaction of the stellar 
wind from a subdwarf star with a planetary co mpanion 



may also cause observable effects ( Geier 2010l) 



10 Conclusions 

We have shown that stellar winds from central stars of 
planetary nebulae resemble stellar winds of luminous O 
stars. In contrast, stellar winds of hot subdwarf stars 
are much weaker and resemble stellar winds of B stars 
in several important aspects. They 



6 



• are driven mainly by lighter elements, 

• show multicomponent effects, 

• may be very hot in outer parts due to inefRcient shock 
cooling. 
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